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Abstract The detection of the orbital period of 4^36 is re- 
ported for the new Intermediate Polar 1 WGAJ1958. 2+3232. 
The orbital period was derived from time-resolved photomet- 
ric and spectral observations. We also confirmed the 733 sec 
spin period of the White Dwarf consistent with the X-ray pul- 
sations and were able to distinguish the beat period in the light 
curve. Strong modulations with orbital period are detected in 
the emission spectral lines from spectral observations. They 
show the presence of a bright hot spot on the edge of the ac- 
cretion disk. The parameters of this recently discovered Inter- 
mediate Polar are determined. 
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1. Introduction. 

Cataclysmic variables (CVs) are close binary systems in which 
mass is transferred from a red dwarf star that fills its Roche lobe 
onto a white dwarf (WD). Intermediate polars (or DQ Her sys- 
tems) are a subclass of magnetic cataclysmic variables with an 
asynchron ously r otating { Psvin < Porb) magnetic white dwarf 



(Patterson, 1994; Warner, 1995). The accretion flow from the 



red dwarf star forms an accretion disk around white dwarf, and 
this disk is disrupted by the magnetic field close to the white 
dwarf. Within the magnetospheric radius, the material is chan- 
neled towards t he ma gnetic polar regions in an arc-shaped ac- 
cretion ( Rosen, |l988l ) 

The recently discovered pulsating X-ray source 
1 WGA J1958.2+3232 (Israel et al., |199S| ) was announced as a 
new Intermediated Polar (IP) by Negueruela et al., ( 2000| ) from 
spectral observations. Strong modulations of this source in 
X-rays were obtained from ROSAT PSPC (721 ± 14 sec) and 
more accurate 734 ± 1 sec from ASCA are presented by Israel 
et al. (1998) and Israel et al. (1999). Photometric observations 



of the optical counterpart of 1 WGA J1958. 2+3232 exhibited 



Table 1. Observations Log 
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strong optical variations, compatible with the X-rays (^ 12 
min) period ( Uslenghi et al., 2000). This modulation was 
interpreted as evidence of a spin period of the WD in this close 
binary system. 

In this paper we present the results of new photometric and 
spectral observations of this system. 



2. Observations. 

The CCD photometric and spectral observations of the 
1 WGA J1958.2+3232 were carried out on 2-5 August 2000 
at the 1.5m and 2.12m telescopes of the Observatorio Astro- 
nomico Nacional, San Pedro Martir of the Institute of Astron- 
omy of UNAM, Mexico. The log of observations is presented 
in Table.|I} 

2.1. Optical photometry 

We obtained i?c-band time-resolved photometry of the optical 
counterpart of 1 WGA J1958.2+3232 during two nights in Au- 
gust 2000 at the 1 .5m telescope. The telescope was equipped 
with a 1024 X 1024 pixel SITE CCD. The frame was reduced in 
size to 450 x 450 pix for faster read-out. It accommodated the 
object and at least two comparison stars in the field of view. The 
exposure times were 120 sec, which leads to the time resolution 
of 169 sec, taking into account dead time between readouts. In 
total the object was monitored during '^13.55h (6.3h the first 
night and 7.25h the second). The data reduction was performed 
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Figure 1. 1 WGAJ1958.2+3232 light curves in Re band ai-e 
presented. Binning time is 120s. 
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Figure 2. The CLEANed power spectrum of Re light curve is 
given. The orbital period ftorb and uJrot spin period frequen- 
cies are marked. In the inset part of the cleaned power spec- 
trum with a fewes number of iterations is shown. The beat fre- 



were corrected for bias and flatfield before aperture photom- 
etry was caiTied out. An estimate of uncertainty of the CCD 
photometry of the optical counterpart of 1 WGA J1958.2H-3232 
was obtained from the dispersion of magnitudes in the differ- 
ential photometry of comparison stars with similar brightness. 
The dispersion ranged from 0.005 to 0.01 mag. We did not ob- 
tain an absolute calibration for our photometric data. 

2.2. Optical spectroscopy 

Time-resolved spectroscopy of the optical counterpart of 
1 WGAJ1958.2H-3232 was obtained on 4-6 Aug, 2000 using 
1 the Boiler & Chivens spectrograph installed in the Cassegrain 
focus of the 2.12m telescope. We used the 4001/mm grating 
with a 13°54 blaze in the second order, combined with the blue 
BG39 filter and CCD TEK1024 x 1024pix with a 0.24/i pixel 
size. The slit width was 1.5 arcsec projected on the sky. This 
combination yielded a spectral resolution of 2.7A FWHM and 
provided a wavelength coverage of A4050 — 5600 Aprom three 
nights of spectral observations the second and third nights were 
disrupted by passing clouds. However the seeing was satisfac- 
tory with images being < 1 .2 arcsec. The slit was oriented with 
position angle of 306° to accommodate nearby star for the flux 
level control. The exposure time in the first two nights was 700 
sec, while on the third night, 700 and 350 sec. The He-Ar com- 
parison spectra were taken every ^120 min. A total of 68 
spectra were obtained. The IRAF long slit spectroscopic reduc- 
tion package was used for extraction of spectra, wavelength and 
flux calibrations. Before head the images were reduced for bias 
and cosmic rays. 

3. Results 

3.1. Light Curve. Orbital and Spin Modulations 

The object demonstrates multi-scale time variability in the 
range 0.3 magnitudes (see. Figjl^). Four pronounced eclipse- 
like depressions obviously shape the light curve. Strong flick- 
ering with optical pulse amplitude (semiamplitude) of about 
0.1 magnitude is also obvious in the light curve detected 
and identified earlier (Israel et al., 1998, Uslenghi et al.. 



quencies corot 
presented. 



n, 



rb, ^rot + ^orb and uj — 1/727.78 sec ai^e 



by using ESO-MIDAS|] and IRAF g softwares. The images 

' ESO-MIDAS is a copyright protected software product of the Eu- 
ropean Southern Observatory, and provides general tools for image 
processing and data reduction. 

^ IRAF is the Image Reduction and Analysis Facility, a general pur- 
pose software system for the reduction and analysis of astronomical 
data. IRAF is written and supported by the IRAF programming group 
at the National Optical Astronomy Observatories (NOAO) in Tucson, 
Arizona. NOAO is operated by the Association of Universities for Re- 



2000) as spin related modulations. The photometric data of 



1 WGAJ1958.2H-3232 was analyzed for periodic ities using 
Discrete Fourier Transform (DFT) code (Deeming, 1975) with 
the CLEAN procedure (Roberts et al., |1987[ ). The CLEANed 
power spectrum (Fig.|]) of Re photometric data shows a clear 
peak at ^orb — 5.54996 ± 0.39624, corresponding to P = 
0.1802 ± 0.0065(i. This peak is caused by above mentioned 
eclipses in the light curve and clearly marks the orbital period 
of the system. 

We found also a significant peak at the spin period ojrot of 
the WD corresponding to 733.82 ± 1.25sec. This period is in 
excellent accordance with that recently discovered by ASCA 
X-ray pulsations (Israel et al., 1999). The beat frequencies at 
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Figures. The diagnostic diagrams for the Hp emission lines. 
The radial velocity semi-amplitude K, the ratio (tk /K are plot- 
ted as a function of the Gaussian separation, obtained for a pe- 
riod of 0.18152 days. 
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Figure 4. The power spectrum of the Re light curve and the 
He II 4686 and H^ radial velocity curves are presented. They 
are scaled to the amplitude of Re power spectrum. The maxi- 
mum peak of frequency corresponds to the orbital period of the 

system Porb = 0.18152 ± O.OOlld . 



l^rot 
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.(,, ojrot + ^orb ^6 also present in the CLEANed 



power spectrum but with a fewer number of iterations (see in- 
sert in the upper right corner of the FigH). The harmonics of 
basic frequencies ilorb and ujrot are detected as well. Besides 
these, there are comparably significant peaks at the periods of 
727.78 sec and 1 .36h. The former was detected also by Uslengi 



et al., (2000) and is probably the one day alias of the uJrot, while 



for the latter we could not find any reasonable explanation. 



3.2. Radial Velocity Variations and Binary System Parameters 

The spectrum of 1 WGA J1958.2H-3232 shows features char- 
acteristic to Cataclysmic Variables. We refer to the Fig.l of 
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Figure 5. The radial velocity curves for the emission lines of 
Hp and He II 4686 phased on the spectroscopic orbital period 
(0.18152 days) are presented upper and middle panels respec- 
tively. The filled circles corresponds to 4 Aug, filled rectangu- 
lars and open deamons are 5 and 6 Aug, respectively. Binning 
time is 700s. Re Ught curve of 1 WGAJ1958.2H-3232 is pre- 
sented in the lower panel. The binning time is 120sec. 




Negueruela et al., (2000), who obtained spectrum of the ob- 



ject in a wider spectral range and with better spectral resolu- 
tion. Meanwhile we obtained time-resolved spectroscopy of 
1 WGA J1958.2H-3232 around the emission lines of H^ And 
He II , covering several orbits. Thus, we were able to consider 
periodical variations in the spectrum of the object, primarily in 
the emission lines. The simple stacking of consecutive spectra 
onto the trailed spectrum showed strong variability in lines. It 
is distinct in Balmer lines and in the higher excitation lines of 
ionized Helium. The Balmer lines are double-peaked with S- 
wave moving inside, which makes it hard to see the periodic 
pattern. In the He II 4686 line the central narrow component 
dominates in most of the phases, and it shows clear sinusoidal 
variation. 

In order to determine the orbital elements we measured the 
radial velocities (RV) of H/3 applying the double Gaussian de- 
convolution method introduced by Schneiger & Young, ( 198C| ) 
further developed by Shafter (1983 ). This method is especially 
efficient for measurements of the orbital motion of CVs with a 
prominent spot at the edge of the accretion disk, contaminating 
the central parts of the emission lines. It allows us to measure 
RV variations using the wings of the lines. The width of Gaus- 
sians were chosen to be slightly larger than our spectral reso- 
lution (8.5 A), where deconvolution was reached at all orbital 
phases. The radial velocities were measured as a function of 
distance a between the Gaussians, and then the diagnostic dia- 
grams were constructed using an initial guess for the orbital pe- 
riod, derived from photometry and from preliminary radial ve- 
locity measurements via Gaussian fits to the lines. The optimal 
value of separation (a — 1175 km /sec ) was determined from 
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Figure 6. The angle of inclination of the system vs. q = 
M2/M1 are shown (solid line). The mass of WD in the system 
Ml vs. q for AI2 = OAlAIgun are given by the dot-dashed line. 
"Best fit" parameter from the secondary star position and the 
ballistic trajectory of the gas stream on the Dopplers maps are 
marked. 

the diagnostic diagrams, and the RV values measured for these 
Gaussian separations were again subjected to a power spectrum 
analysis in order to refine the period. The spectroscopic period 
peaked at slightly longer value, than the photometric (however 
within errors of the photometric period). This method quickly 
converged and after two iterations no further improvement was 
achieved. The diagnostic diagrams for H^ are shown in Fig.^. 
Fig.||(top) shows the H^ radial velocity curve. 

A narrow single Gaussian profile was fitted to the promi- 
nent emission features in the profile of He II A4686A, and the 
measured line centers were used to determine the radial ve- 
locity solution for 1 WGA J1958.2H-3232. The radial velocity 
curve for He II A4686 line is presented on the middle panel of 
Fig.|l These measurements were also subjected to the power 
spectrum analysis. The obtained orbital period is in good agree- 
ment with the values derived from the photometry and the 
H^ line radial velocities. The power spectra around the val- 
ues corresponding to the orbital period from these three in- 
dependent determinations are plotted in the Fig.0. One can 
see the excellent match of the central peak. We adopted the 
0^18152 ± 0'?00011 as the final value for the orbital period 
of 1 WGAJ1958.2H-3232 from our observations. Longer time 
base observations are needed to improve this value. 

Each of the radial velocity curves was fitted using a least- 
squares routine of the form 



v{t) =-fo + Ki* sin(27r(i - to)/P), 



(1) 



where -fo is the systematic velocity of the system, and Ki is 
the semi-amplitude of radial velocity, both in km s^^. The ob- 
servation time is t, the epoch to = 2451761^2281 ± 0^0001 
corresponds to the +/- zero crossing of the Hp radial veloc- 
ity curve, and therefore is a superior conjunction of the binary 



system (secondary located between observer and the WD). Ac- 
cordingly the phase 0.0 was calculated and operated further at 
this epoch. Table.Hgives summary of the radial velocity fits for 
H^ and He II 4686 emission line. 

Table 2. Radial velocity solution parameters 



Name 



7o 
km/s 



km/s 



HeIM686 



-39 ±5 
-70.2 ±15.0 



74 ±7 
197.4 ±25.7 



After refining the orbital period from spectroscopy, and de- 
termining the phase 0.0, the photometric light curve was folded 
by the corresponding parameters and presented in the lower 
panel of Fig. 1} 

Several conclusions can be made after considering the three 
curves in Figjs] in conjunction with and taking into account 
common knowledge of Intermediate Polar systems (Patterson, 



1994, Warner, 19951 



- The double peaked Balmer lines are due to the presence of 
an accretion disk or ring orbiting the primary WD in this 
IP. RV variations in the wings of lines describe the rotation 
of the primary of the binary system. 

- The S-wave present in the Balmer and He II emission lines 
is evidence for a hot compact region on the accretion disk. 

- From the difference of phases and amplitudes of radial ve- 
locities of wings of H/3 and the naiTow component of He II 
we can localize the hot spot in the second quarter of the 
disk, if the center of coordinates is at the WD. It is the 
usual location of the spot originating from the impact of 
the mass transfer stream with the accretion disk. However 
it is also the area which most commonly heated by the ener- 
getic X-ray beam from the magnetically accreting pole on 
the surface of the WD in IPs (see the sketch in Hellier et 
al.,p989|). 



- The hot spot is self eclipsed by the accretion ring as follows 
from the phasing of the light curve, the dips in the light 
curve centered at the phase « 0.12cf)orb- The primary at this 
phase starts to move toward the observer and the hot spot 
is on the opposite side of the ring approaching maximum 
velocity. 

We estimate the mass and radius of the secondary star as 
M2 = O.41M0and R2 = 0A7Rq from 

M2/Msun = 0.0751P(/l)l■l^ 

R2/Rsun = 0.101P(/^)l■"^ 1.4 < P{h) < 12 (2) 

the mass-period and radius-period relations of EchevaiTia, 



(1983) 



The mean mass estimate of 76 CV white dwarfs is Mwd = 
O.86M0 (Sion, |1999[ ). Webbink, ( |1990| ) gives statistically av- 



erage white dwarf masses ratios (q = 0.29) and average masses 
for all systems (Mwd = 0.61Mq) below the period gap and 
(q = OMMivd = 0.82Mq) above the period gap. 
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Figure 7. Trailed, continuum-subtracted, spectrums of IWGA J 1958.2+3232 plotted in two cycles. Doppler maps of the emission 
lines He II (left panel), Hp (right panel) in velocity space (Vx, Vy) are given. A schematic overlay marks the Roche lobe of the 
secondary, the balUstic trajectory and the magnetically funneled horizontal part of the accretion stream. The secondary star and 
gas-stream trajectory are plotted for K=74 km/s and q=0.46. 



From our spectroscopic radial velocity solution, we can de- 
termine preliminary values for the basic system parameters of 
1 WGAJ1958. 2+3232. The binary orbital plane inclination can 
be determined from (see Downes et al., 1986; Dobrzycka & 
Howel,|l992l|: 



. 3,, KfP .q 



1 



2itGM2 ' q 



(3) 



if the mass ratio of the system is known. The dependance of 
i versus q — M2/M1 in the range 0.4 up to 0.75 is shown in 
Fig.g for the above determinated Ki h„ ■ 

We attempted to refine this estimate for 
1 WGAJ1958. 2+3232 by constraining Doppler tomograms 



Table 3. The 1WGA1958. 3+3232 adopted system parameters 



Parameter 


Value 


Parameter 


Value 


Porb 


0.18152d 


R2 


0.47 Rq 


Prot 


733.7s 


a 


1.5 Rq 


M2 


0.41AfQ 


i 


35° 


q 


0.46 


Afi 


0.9Mq 



from observed emission line profiles. Doppler tomography 
is a useful tool to extract further information on CVs from 
trailed spe ctra. T his method, which was developed by Marsh 
& Home, (1988 1, uses the velocity profiles of emission-lines 
at each phase to create a two-dimensional intensity image in 
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Figure 8. Doppler maps of the emission line blend C Ill/N III 
in velocity space (14, Vy) is given. A schematic overlay marks 
the Roche lobe of the secondary, the ballistic trajectory and the 
magnetically funneled horizontal part of the accretion stream. 
The secondary star and gas-stream trajectory are plotted for 
K=74 km/s and q=0.46. 

velocity-space coordinates (Vx,Vy). Therefore, the Doppler 
tomogram can be interpreted as a projection of emitting 
regions in cataclysmic variables onto the plane perpendicular 
to the observers view. We used the code developed by Spruit 



(1998) to constrain Doppler maps of 1 WGAJ1958.2H-3232 



with maximum entropy method. The resulting Doppler maps 
(or tomogram) of emission lines of H/3 , He II and the blend 
of C Ill/N III are displayed as a gray-scale image in Fig.0 
and Fig.^. Also in Fig.0 are displayed trailed spectra of H/? 
and Hell in phase space and their corresponding recon- 
structed counterparts. Two features in the maps are distinct: 
an accretion disk seen as a dark circle extending to up to 
^ —700 km/sec on H/3 doppler tomogram and a bright spots 
detected in all three maps to the left and below the center of 



mass at velocities T4 ^ —225 km s , Vy 



-100 km s in 



He II 4686. Apart from the spot a cometary tail linked to it and 
extending to T4 w Vy « —500 km s^^ can be clearly seen on 
the He II map. These we identify with the mass transfer stream 
and its shape was essential for our selection of the ballistic 
trajectory. A helpful assistance in interpreting Doppler maps 
are additional inserted plots which mark the position of the 
secondary star and the ballistic trajectory of the gas stream. 
Here we used our estimates of Porb and M2 with various 
combination of i and q from Fig.H in order to obtain the "best 
fit" (by simple eye inspection) of calculated stream trajectory 
with the gray scale image. Our preferred solution of inclination 
is i = 35° deg. It is marked on Fig.pl and is given in TablcB. 
Of course other close solutions are applicable. Comparing the 
location of spots in Doppler maps of He II and Hf3 we can 



distinguish actually two hot spots on the disk (Fig.^. The 
elongated spot coinciding in both emission lines is caused 
probably by the mass transfer stream and the shock of impact 
with the disk, while the compact dense spot toward negative 
Vy's, much better seen in H/3 than in the two other lines, is 
a result of heating of the disk by X-ray beam. The C Ill/N III 
pattern mostly repeats that of He II with less intensity of course 
(Fig.|). 

4. Conclusions 

The 1 WGA J1958.2H-3232 results to be a classical "textbook" 
Intermediate Polar. It has an orbital period Porb = 4^36 above 
the period gap, as the vast majority of IPs. It exhibits X-ray 
and optical coherent pulsations of the order of ~ 0.05Por&, 
undoubtfully originating from asynchronous spin of magnetic 
WD in a close binary system. Also the beat period in optical 
light is detectable. This is another characteristic of Intermedi- 
ate Polars. 

Other orbital parameters derived from assumption that the 
system obeys the Por6^M2 relation for CVs are also in-line 
with acc umula ted data on oth er IPs and theoretical aspects 
(Warner, |l995i Paterson, |l994 see also URL0 ). 

The radial velocity curves, the light curve and the Doppler 
tomography confirm the presence of an accretion ring around 
the WD and the existence of hot spots caused by heating of 
parts of the disk by the X-ray beam and from interaction with 
the mass transfer stream. 
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